Instability of condensation in the zero-range process with random interaction 
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The zero-range process is a stochastic interacting particle system that is known to exhibit a 
condensation transition. We present a detailed analysis of this transition in the presence of quenched 
disorder in the particle interactions. Using rigorous probabilistic arguments we show that disorder 
changes the critical exponent in the interaction strength below which a condensation transition may 
occur. The local critical densities may exhibit large fluctuations and their distribution shows an 
interesting crossover from exponential to algebraic behaviour. 
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The zero-range process is a stochastic lattice gas where 
the particles hop randomly with an on-site interaction 
that makes the jump rate dependent only on the local 
particle number. It was introduced in [T| as a mathemat- 
ical model for interacting diffusing particles, and since 
then has been applied in a large variety of contexts, of- 
ten under different names, (see e.g. [2J and references 
therein) . The model is simple enough for the steady state 
to factorize, on the other hand it exhibits an interesting 
condensation transition under certain conditions. Viz. 
when the particle density exceeds a critical value p c the 
system phase separates into a homogeneous background 
with density p c and all the excess mass concentrates on 
a single lattice site. This has been observed and studied 
in some detail in experiments on shaken granular media 
d, [I], and there is a well-established analogy to Bose- 
Einstein condensation @, H, Q . It is also relevant as a 
generic mechanism for phase separation in single-file dif- 
fusion Q and condensation phenomena in many complex 
systems such as network rewiring or traffic flow @ , for 
a review see Q • 

The transition can be caused by site-dependent jump 
rates g x [f| due to the slowest site acting as a trap. It 
also appears in a more subtle fashion in homogeneous 
systems where condensation may result from the growth 
of large clusters on the expense of small clusters, if the 
jump rates g{n) as a function of the number n of particles 
on the starting site have a decreasing tail. Such a model 
with a generic power law decay 



j(n) = 1 + b/n a 



(1) 



with positive interaction parameters b, a has been intro- 
duced in [T(|. Condensation occurs if < a < 1 and 
b > or if a = 1 and b > 2. The condensation transition 
is understood also rigorously in the context of the equiv- 
alence of ensembles (Til [t3 |. and more recently many 
variants of Q} have been studied @, @, [U, [T3, GB El] . 
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FIG. 1: (Color online) Change of the phase diagram under 
random perturbations of the jump rates ([T|l. Disorder changes 
the critical interaction exponent from 1 to 1/2 and leads to 
a critical density that depends on the system size. Inside 
the red-shaded region condensation occurs above a non-zero 
critical density for b > 0, and for b < (see ([2])) and negative 
a the critical density is zero. 



All previous studies of zero-range processes assume 
that the interaction between particles is strictly equal 
on all sites. In this paper we study the effect of dis- 
order on this interaction and show that even a small 
random perturbation of the n-dependence of the jump 
rates g(n) leads to a drastic change in the critical be- 
haviour. Namely, for positive b condensation only occurs 
for < cr < 1/2 (see Figure [1]), excluding in particular 
the frequently studied case a = 1. Moreover, the critical 
densities are site-dependent random variables with non- 
trivial distributions and depending on the parameter val- 
ues they may exhibit large fluctuations. Due to the wide 
applicability of the zero-range process, the change of the 
critical interaction exponent a is particularly relevant for 
applications as is explained in the conclusion. 

We consider a lattice Al, which we take to be periodic 
and of finite size |AjJ = L. A configuration is denoted by 
(j]x)xeA where r/ x G {0, 1, . . .} is the occupation number 
at site x. The dynamics is defined in continuous time, 
such that with rate gxi'Hx) site x G A^ loses a particle, 



2 



which moves to a randomly chosen target site y according 
to some probability distribution p(y—x). For example 
in one dimension with nearest neighbour hopping, the 
particle moves to the right (left) with probability p (1— p). 

A generic perturbation of the jump rates |T]) can be 
additive or multiplicative, but since the condensation be- 
haviour is determined only by the tail of the jump rates 
for large n, both choices are essentially equivalent. They 
can be written in a convenient general way, 

9x (n) = e E ^ for n > 1 , g(0) = , (2) 

where the exponents are given by 

E x {n) = e x (n) + b/n <J , bet, a > , (3) 

with e x (n) being iid random variables with respect to x 
and n. Without ^-dependence this would amount merely 
to a change of the function g which might be interesting, 
but is a degenerate problem in terms of generic perturba- 
tions. On the other hand, the effect of spatially inhomo- 
geneous jump rates favouring condensation on slow sites 
has already been studied [f|. Therefore we concentrate 
on disorder with spatially uniform mean M(e x (ri)) = 
and variance S 2 > 0, in order to focus on the basic nov- 
elty which is the supressing effect on condensation for a 
generic perturbation of (TT]) . For the same reason we have 
chosen the jump probability p(y—x) to be homogeneous, 
since spatial dependence there leads to the same effect as 
spatially inhomogeneous jump rates p7| . 
Note that for e x (n) = (i.e. 6 = 0) the asymptotic be- 
haviour of g x {n) is given by ([1]). All data shown in this 
paper are for uniform e x (n) ~ U(—a, a), characterized 
by S 2 = a 2 /3. But our analytical results are of course in- 
dependent of the distribution of the perturbation as well 
as the exact form of the jump rates ([2]). For negative a 
the rates (J2J) are increasing in n for positive b and hence 
there is no condensation. For negative b the rates tend to 
zero, which means that there is condensation at critical 
density p c = 0. This is an essentially trivial feature of the 
model which is robust against perturbation by disorder. 
We therefore focus on positive interaction exponent a. 

It is well known (see e.g. Q, [TH) that the process has 
a grand-canonical factorized steady state with single- 
site marginal 

vxA n ) = — rr ll^^r 1 = 

1 ™ 

= ^) exp ( n/i -g^ (fc) ) ' (4) 

where the chemical potential fi £ R fixes the particle 
density. This holds independently of the distribution of 
target sites p(y — x) and for each realization of the e x (fc), 
i.e. i/p is a quenched distribution. The single-site nor- 
malization is given by the partition function 

oo n 

z x{p) = ^2cxp(?in-^2E x (k)) (5) 

n=0 k=l 
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FIG. 2: (Color online) Typical realizations of y/n£ x (n) 
(drawn with 8 2 — l) leave the area enclosed by the dashed 
parabola only finitely many times, but cross the full line 
-j^ n 1 ~ CT for a > 1/2 infinitely often (here a = 0.75, b = 1). 

which is strictly increasing and convex in /i The 
local density can be calculated as usual as a derivative of 
the free energy 

P W = {Vx) = (6) 

and is a strictly increasing function of /i. By (..) we 
denote the (quenched) expected value with respect to 
for fixed disorder, i.e. fixed realization of the e x (k). To 
study the condensation transition we have to identify the 
maximal or critical chemical potential fi c £ R, such that 
z x (fi) < oo for all [i < [i c . 

The single-site marg inal (gl) is a function of the disorder 
and for b = it has the distribution of a geometric ran- 
dom walk with deterministic drift fi. The critical chemi- 
cal potential in this case is simply /i c = and as n — > oo 
the distribution converges to a log- normal. With 

non-zero b the drift term becomes n-dependent and a 
more detailed analysis is required. With Q we have to 
leading order as n — > oo 

fc=i k 

where by the central limit theorem 
1 n 

Un) — T-foY; e *( fc ) ™ ^ (°> !) ( g ) 
v n k=i 

converges to a standard Gaussian. Moreover, the process 
(\/n£, x (n) :ij£N) is a random walk with increments of 
mean zero and variance 1. Since the fluctuations of such 
a process are of order ^/n we have for all C £ R 

F(t; x (n) ^ C for infinitely many n) = 1 , (9) 

and for all 7 > 0, C > 

P(|^j;(n)| > Cn 1 for infinitely many n) = . (10) 

This is a direct consequence of the law of the iterated 
logarithm (see e.g. [lj|, Corollary 14.8) and is illustrated 
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in Figure [2j Together with ([5]) and (JT]) this implies that 
z x (p) < oo for all /i < with probability one. So for 
almost all (in a probabilistic sense) realizations of the 
e x (k) the critical chemical potential is \x c = and 



OO / Tl 

Zx(p-c) = E cxp ( ~ E ( e *( fc ) 



n=0 



fc=l 



6 



(11) 



For certain values of a and 6, z x (p, c ) < oo is possible and 
v x,ii a can be normalized, which is a necessary condition 
for a condensation transition [13] • We find 

b < 0: In this case (JJJ) and (J9j) imply that there are 
infinitely many terms in which are bounded below 
by 1. Since all terms of the sum are non-negative it 
diverges and z x {p c ) = oo with probability one. 

b > 0: In this case the asymptotic behaviour of the 
terms in the sum (fTTj) depends on the value of a > 0, 
since the sign of the exponent can change. 

• For a > 1/2, n 1 ~ a <C s/n and Q is dominated 
by 5^/nt] x {n). Applying (JOj) with C = we get 
z x {pc) = oo with probability one. 

• For (7=1/2 both terms in ([7]) are of the same order 
~ Jn and 



since n 



fe=i 



-<V^(U«)+t) 



(12) 



Again, j9} this time with C = 2b/ 6 implies 
ZxifJ'c) = oo with probability one. 



• For < cr < 1/2 we have n 1 ^ a > y/n and is 
dominated by j^-n 1- " 7 . We apply iflQ]) for 7 = 
1 — cr — 1/2 > Oto see that the random quantity 
£ x (n) can change the sign of the exponent in (TTT]) 
only for finitely many terms in the sum. Therefore 
z x (Pc) < 00 with probability one since jzr^n 1 ~ a 
has a fixed negative sign in ((TTJ) . 

Whenever z x (p c ) = 00 the local critical density 

pl~p x {n c ) = oo (13) 

also diverges and there is no condensation transition (see 
e.g. [ll|, Lemma 1.3.3). But for b > and < a < 1/2 
we have z x (^, c ) < 00 and by the same argument as above 
it follows that 



Pc 



z x (n, 



- OO 



-£L=l(e*(fc)+7^) 



< OO 



(14) 



with probability one, since the factor n in the sum only 
gives a logarithmic correction in the exponent. Therefore 
there is condensation since the grand-canonical product 
measures only exist up to a total density of 



Pc 



(L) := ~ £ Pi , 



(15) 



xeA L 



which depends on the e x (n) and the size of the lattice L. 
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FIG. 3: (Color online) Background density pb g (|17|) as a 
function of N/L for b = 1.2. Data for fixed disorder with 
8 2 — 1/12 and L = 1024 show condensation for a — 0.2 and 
no condensation for a = 0.8. Data for S 2 — are shown by 
unfilled symbols, dotted lines indicate the critical densities. 
Data points are averages over 100 MC samples with errors of 
the size of the symbols. 



If the actual number of particles N is larger than 
L Pc{L), all sites except the 'slowest' one contain on aver- 
age (t] x )n = Pc particles. They form the so-called critical 
background, since their distribution has non-exponential 
tails. By (..)jv we denote the (canonical) expectation 
conditioned on the total particel number TV. The slowest 
site, say y, is defined by p v c > p x c for all x 7^ y. By the 
conservation law it is required that 



(Vy) 



N 



N ■ 



J2p x c = o(n) , 

x=£y 



(16) 



i.e. it contains of order N particles and forms the conden- 
sate. This interpretation is in accordance with previous 
results and has been proved rigorously in poj in the limit 
N —> 00 for the unperturbed model. For the perturbed 
model we support this conclusion by MC simulations, 
some of which arc shown in Figure [3] For fixed L we plot 
the stationary background density 



Pbg 



y)Nj 



(17) 



as a function of the total density p = N/L. For a = 0.2, 
Pbg converges to a critical density p c {L) which is slightly 
higher than for the unperturbed model. The overshoot 
of pbg for densities close to p c (L) is due to sampling from 
the canonical rather than the grand canonical ensemble, 
which has been observed already in (22j. For a = 0.8, pt, g 
increases approximately linearly with p, which is clearly 
different from the unperturbed model condensing with 
critical density 2.77. 

For the perturbed system the critical density p c (L) is 
a random variable, which according to (fTB"|) converges in 
the thermodynamic limit to the expected value E(Pc) 
with respect to the disorder. This is in general hard to 
calculate [23|, even for simple choices of the e x (n). De- 
tailed numerical estimates of the distribution of p x c (some 
of which are shown in Fig. 0|) indicate that depending 
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FIG. 4: (Color online) Cumulative tail of the distribution of 
from 5- 10 5 independent numerical calculations of (|14p with 
6 = 1, where infinite sums have been cut off at n = 10 6 . The 
full line corresponds to a cumulative exponent —1. The data 
exhibit power law decays and for some parameters also large 
exponential cut-offs, e.g. ~ e~ r//1 ° for a = 0.45, 8 2 = 1 fitted 
by a dashed line. The expectations are finite but may be very 
large, E(p?) = 1.60 (A), 2.55(f), 1106 (x). 

on the system parameters the cumulative tail is either 
algebraic or consists of an algebraic part with an expo- 
nential cut-off at large values. This can be explained 
heuristically by the interplay of the two terms in j7]) that 
determine the main contributions to the partition func- 
tion [24|. The exponents of the purely algebraic tails is 
smaller than but often close to —1 and the length scale 



of the exponential tails can be extremely large (cf. Fig. 
|4|. Therefore p x has a finite mean which determines the 
thermodynamic limit of (fT5|) 

Pc := Urn p c (L) = E(p x c ) < oo , (18) 

L — >oo 

but for finite L the p x can exhibit large fluctuations, re- 
sulting in high values for p c and slow convergence of (|18|1 . 

The interaction encoded in the jump rates of a zero- 
range model represents an effective interaction for which 
space-dependence and randomness due to microscopic 
impurities or heterogeneities in the case of complex sys- 
tems has to be taken into account. We have shown that 
generic interaction disorder reduces the critical interac- 
tion exponent from a = 1 to a = 1/2. In particular, 
this implies that for the most-studied case a = 1 there is 
no condensation transition in the presence of interaction 
disorder. This case is relevant for the mapping of the 
zero-range process to exclusion models where it becomes 
an effective model for domain wall dynamics and there- 
fore a powerful criterion for phase separation in more 
general particle systems Q ■ In this mapping the spatial 
interaction disorder maps into hopping rates which de- 
pend in a random fashion on the interparticle distance. 
Then our results imply that in such heterogeneous finite 
systems the change of the critical interaction exponent 
has to be taken into account and finite size effects may 
play a major role due to large fluctuations of local critical 
densities. 
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